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Abstract 

Eukaryotes have evolved complex defense pathways to combat invading pathogens. Here, we investigated the role of the 
Arabidopsis thaliana heterogeneous nuclear ribonucleoprotein (hnRNP-Q) LIF2 in the plant innate immune response. We 
show that UF2 loss-of-function in A. thaliana leads to changes in the basal expression of the salicylic acid (SA)- and jasmonic 
acid (JA)- dependent defense marker genes PR! and PDF1.2, respectively. Whereas the expression of genes involved in SA 
and JA biosynthesis and signaling was also affected in the Iif2-1 mutant, no change in SA and JA hormonal contents was 
detected. In addition, the composition of glucosinolates, a class of defense-related secondary metabolites, was altered in 
the Iif2-1 mutant in the absence of pathogen challenge. The Iif2-1 mutant exhibited reduced susceptibility to the hemi- 
biotrophic pathogen Pseudomonas syringae and the necrotrophic ascomycete Botrytis cinerea. Furthermore, the Iif2-1 sid2-2 
double mutant was less susceptible than the wild type to P. syringae infection, suggesting that the Iif2 response to 
pathogens was independent of SA accumulation. Together, our data suggest that Iif2-1 exhibits a basal primed defense 
state, resulting from complex deregulation of gene expression, which leads to increased resistance to pathogens with 
various infection strategies. Therefore, LIF2 may function as a suppressor of cell-autonomous immunity. Similar to its human 
homolog, NSAP1/SYNCRIP, a frans-acting factor involved in both cellular processes and the viral life cycle, LIF2 may regulate 
the conflicting aspects of development and defense programs, suggesting that a conserved evolutionary trade-off between 
growth and defense pathways exists in eukaryotes. 
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Introduction 

Living orgaiiLsms are exposed to various abiotic and biotic 
stresses and continuously integrate diverse signals to limit tlie 
spread of microbial pathogens. RNA regulatory processes, which 
modulate transcriptional programs and post-transcriptional events 
in response to various cues, are an important component of the 
tolerance and adaptation strategies of organisms. Thus, RJMA- 
binding proteins (RBPs), which are involved in various aspects of 
RNA processing (i.e., mRNA maturation, editing, splicing, and 
mRNA trafficking), are thought to be key regulators of stress 
responses in both animals and plants. However, this possibility has 
been scarcely studied. 

The heterogeneous nuclear ribonucleoprotein (hnRNP) group 
of RBPs forms a large family of ancestral modular proteins with a 
high degree of functional diversification [1]. In humans, this group 
contains 37 members, only a few of which are involved in the stress 
response. For instance, hnRNP-K participates in the response to 
genotoxic stress [2]. The human NSl -associated protein 1/ 
synaptotagmin-binding cytoplasmic RNA-interacting protein 



(NSAPl/SYNCRIP), another ceUular hnRNP protein, partici- 
pates in the translational activation of both cellular and viral 
RNAs and functions in the Hepatitis C virus (HCV) life cycle by 
interacting with HCV RNA [3,4,5]. Knock-down of NSAPl/ 
SYMCRIP significandy decreased the amount of HCV RNA in 
mammalian cells, suggesting that NSAPl/SYNCRIP is a negative 
regulator of viral defense responses. Recently, hnRNP-I was 
shown to interact with the long non-coding RNA RoR to modulate 
the expression of the tumor suppressor p53 in response to DNA 
damage and to participate in a surveillance network [6]. 

Plants contain more RBPs than do animals, but the number of 
plant hnRNPs is similar in both plants and animals. The functions 
of hnRNPs remain poorly described in plants [1,7]. Few studies 
have reported a role for plant hnRNPs in the stress response [8]. 
Until now, the glycine-rich RBP AtGRP7 was the only hnRNP 
known to play a role in the plant's response to pathogens [9]. 
AtGRP7 is the substrate of the type III effector HopUl, which is 
injected by the bacterial pathogen Pseudomonas syringae into plant 
host cells. HopUl mono-ADP ribosylates a conserved arginine 
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residue of AtGRP7, preventing binding to immunity-related RNA 
[10,11]. 

In a previous study, we identified LIF2, a novel hnRNP of the 
cruciferous plant Arabidopsis thaliana [12]. LIF2 is a nucleo- 
cytoplasmic RBP that contains three RNA Recognition Motifs 
(RRMs), the most frequent RNA-binding domain in the hnRNP 
family. LIF2 and its close homologs, AtLILl and AtLIL2, belong 
to the hnRNP-Q subfamily and LIF2 is structurally homologous to 
human NSAP1/SYNCRIP[12]. LIF2 interacts in vivo with LHPl, 
aPolycomb Repressive Complexl (PRCl) subunit [13,14]. LIF2 is 
involved in the maintenance of plant cell identity and cell fate 
decision [12]. Indeed, loss-of-function of LIF2 affected various 
aspects of growth and development, such as flowering time and 
leaf size. More dramatically, lif2 mutations induced indeterminate 
growth of ovaries, resulting in the formation of ectopic inflores- 
cences bearing severely affected flowers. Besides its developmental 
functions, LIF2 might be involved in the stress response. Indeed, a 
large set of stress-related genes (193) was found to be deregulated 
in li/2-1. The set of common genes deregulated in both lif2-l and 
Ihpl mutants was even more enriched in stress-related genes than 
was that in liJ2-l [12]. Furthermore, loss-of-function of LHPl (also 
named TFL2 or TU8) altered the glucosinolate profile and reduced 
symptoms in response to infection by the obligate biotrophic 
fungus Plasmodiophora brassicae, which causes clubroot disease, a 
damaging disease in Brassicaceae [15,16]. 

These data suggest that LIF2 and LHPl have common 
functions in the stress response and prompted us to investigate 
the function of LIF2 in biotic stress responses. Here, \vc show that 
lif2 mutations conft^rred altered susceptibility to the necTotrophic 
fungal pathogens Sclerotinia sclerotiomm (S. sclerotiorum) and Bohytis 
cinma {B. cinered) and to the hemi-biotrophic bacterial pathogen P. 
syringae pv. tomato [Pst), three biological agents with substantial 
impacts on the agronomical production of various plant species 
[17,18,19]. To better understand the altered pathogen suscepti- 
bility of the lif2 mutants, we investigated the expression of defense 
marker genes regulated by jasmonic acid (JA) and salicylic acid 
(SA), two critical signaling hormones in the activation of plant 
defense, as well as the production of stress-associated metabolites 
and hormones. Finally, we showed that the transcriptomic profiles 
of stress-response factors in liJ2 are specific and we identified key 
defense regulators (such as the WRKY18 and WRKY33 transcrip- 
tion factors, which are known to be involved in B. cinerea 
susceptibility), whose deregulation may contribute to the observed 
pathogen response. Together, our data suggest that the LIF2 
hnRNP-Q^ may suppress the plant immune response by an 
unknown SA-independent pathway. Given that the human 
homolog of LIF2, the NSAPl/SYNCRIP protein, is a frmj-acting 
factor involved in both cellular processes and the viral life cycle 
[3,4,5], we propose that the conserved hnRNP-Q_ proteins may 
have an evolutionary conserved function in regulating the trade-off 
between growth and defense in eukaryotes. 

Results 

Analysis of the Iif2-1 transcriptome reveals a potential 
stress-related function for LIF2 

In a previous transcriptome profiling experiment, we showed 
that genes that were deregulated in the A. thaliana Iif2-J nuU 
mutant were greatly enriched in Gene Ontology (GO) terms 
involved in responses to stress stimuli [12]. Here, we further 
examined the deregulated gene set (1008 genes) using Singular 
Enrichment Analysis (SEA) implemented in the agriGO toolkit 
[20]. This analysis revealed 293 deregulated genes associated with 
the GO term "response to stimulus" and 193 with "response to 



stress" (Table 1). We noticed that the normixi fri;quency (NF) of 
the "response to biotic stimulus" (GO:000960, NF 4.22, p-value 
7xlO~^^) was higher than that of the "response to abiotic 
stimulus" (GO:0009628, NF 2.75, p-value 2x10""^) (Table 1). 
Furthermore, several genes were associated with GO terms related 
to the JA-defense signaling pathway, the glucosinolate metabolic 
pathway, and responses to fungal and bacterial pathogens 
(Table 1). These data suggest that even if LIF2 is induced only 
weakly upon pathogen infection (Fig. SI), this gene might be 
involved in defense responses to various bioaggressors. 

JA and SA pathways are altered in the Iif2-1 mutant 

Based on the lif2 transcriptome data [12], we further analyzed 
the expression of a key marker gene of the JA-mediated defense, 
PDF1.2 (a JA-responsive gene), and the expression of genes of the 
JA biosynthesis and signaling pathways, L0X3, AOS, AOC, and 
0PR3 (genes encoding enzymes of the JA biosynthesis pathway), 
JARl (encoding an enzyme that converts JA to the bioactive JA-Ile 
molecule), and COIl (encoding an F-box subunit of the JA- 
receptor complex) [21] (Fig. lA). All of these genes were 
downregulated in lif2-l, suggesting that the JA defense pathway 
is globally repressed in the lif2-l mutant. 

Since the JA-dependent defense network cross-communicates 
with the SA-dependent signaling pathway to fine-tune the plant's 
defense response [22,23], we assessed the expression of ICSl/ 
SID2, which is involved in SA biosyntiiesis; ENHANCED DISEASE 
SUSCEPTIBILTiri {EDSl) and PAD4, which are involved in die 
SA-signaling pathway; and NONEXPRESSOR OF PR GENES 1 
(JVPRl) and PATHOGENESIS-RELATED 1 {PRl), which are 
important positive regulators of SA responses [24,25]. EDSl, 
PAD4, and ICSl were upregulated in the lif2-l mutant relative to 
wild-type plants (Fig. IB). PRl was also upregulated in lif2-l, 
whereas MPRl expression was not significantly affected in lif2-l 
(Fig. IB). These data suggest that most genes involved in the SA- 
mediated defense pathway are constitutively upregulated in the lif2 
mutants, in the absence of pathogen. The upregulation of genes of 
the SA signaling pathway and the downregulation of genes of the 
JA signaling pathway in lil2 illustrates the antagonistic action of 
these two pathways. 

Prompted by the finding that genes involved in JA- and SA- 
defense-related networks are differentially expressed in lif2-l, we 
quantified the endogenous levels of these hormones in lif2-l using 
HPLC-electrospray-MS/MS (Fig. IC). Surprisingly, whereas ICSl 
was upregulated in lif2-l, free SA levels were similar in the mature 
rosette leaves oi lif2-l and wild-type plants (Fig. IC). Similarly, the 
JA levels were not significantly affected in the liJ2-l mutant 
(Fig. IC). Furthermore, we found that the le\'els of auxin (AIA) and 
abscissic acid (ABA), two other phytohormones recently shown to 
be involved in the stress response [26,27], were not significantiy 
altered in lif2-l in the absence of pathogen challenge (Fig. IC). 

The Iif2 mutants are less susceptible to the bacterial 
pathogen Pseudomonas syringae 

As JA- and SA-defense-related gene expression were altered in 
the lif2 mutant (Fig. lA-B), we investigated the response of lif2 to 
inoculation with the hemi-biotrophic bacterial pathogen P. syringae 
[28]. The bacterial growth of the virulent DC3000 strain was 10 
times lower in lif2-l leaves than in wild-type leaves at 24 hours 
post-inoculation (hpi) (Fig. 2A). Consistently, the disease symptoms 
were reduced on the leaves of two independent lines, lif2-l and 
lif2-3, whereas severe symptoms were observed on the leaves of 
wild-type plants five days post-inoculation (dpi) (Fig. 2B). As 
expected, the lif2-c line, a lif2-l mutant line complemented with 
LIF2 under the control of its own regulatory regions, was as 
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Figure 1. JA- and SA-dependent signaling pathways in the Iif2-1 mutant. (A) Expression levels of genes Involved in the JA-dependent 
signaling pathway {L0X3, AOS, AOC, 0PR3, COIl, and JAR1) and of a JA-responsive marker gene {PDF1.2) in the rosette leaves of seven-week-old wild- 
type (WT) and Iif2-1 plants. (B) Expression levels of genes of the SA-dependent signaling pathway (EDSl, PAD4, and ICS1) and of SA-responsive marker 
genes (NPR1 and PRl) in seven-week-old wild-type (WT) and Iif2-1 rosette leaves. CC-NB-LRR and TIR-NB-LRR are disease resistance proteins with 
coiled-coil (CC), nucleotide-binding (NB), leucine-rich repeat (LRR), or Toll-lnterleukin Receptor (TIR) domains. (C) Quantification of phytohormone 
contents in the rosette leaves of three-week-old plants using HPLC-electrospray-IVlS/IVlS. The amount of JA was expressed as a ratio of peak areas 
(209>62/214>62) per fresh weight (FW). The amount of other hormones was expressed in ng/FW. The bars represent standard deviation. 
doi:1 0.1 371/journal.pone.0099343.g001 



susceptible as the wild-type plant (Fig. 2B). Similar results were 
obtained with the avirulent DC3000 avrRpml strain (Fig. 2C). 
These data are consistent with the observation that SA-related 
genes are constitutively expressed in the lif2 mutant. 

In response to P. syringae infection, the JA-related biosynthesis 
genes were upregulated in wild-type plants and to a lesser extent in 
lif2-l (Fig. 3), with the exception oiJARl, which was not induced. 
coil, which encodes a component of the JA-receptor complex 



interacting with the bacterial phyto toxin coronatine, was upregu- 
lated in both genotypes (Fig. 3A). PDF1.2 was downregulated in 
wild-type plants and not activated upon infection in lif2. In 
contrast, genes involved in SA biosynthesis and signaling pathways 
were upregulated in both wild-type and lif2 plants upon DC3000 
inoculation (Fig. 3B). PAD4 induction was weaker in lif2 than in 
wild-type plants, whereas PRl expression was similar in the 
mutant and wild-type plants upon P. syringae infection. Therefore, 
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Figure 2. The ///? mutants are less susceptible to P. syringae \niect\on. (A) Bacterial growth of the virulent DC3000 strain in wild-type (WT) and 
Iif2-1 rosette leaves at 24 hours post-inoculation (hpi). (B) Rosette leaves imaged 5 days post-inoculation (dpi) with the virulent DC3000 strain. The 
Iif2-1 and lif2-3 plants had similar responses, whereas the complemented Iif2-1 mutant (//72-c) behaved similarly to WT plants. Four independent 
experiments were performed with similar results. (C) Bacterial growth of the avirulent DC3000 avrRpml strain in rosette leaves at 24 hpi. For bacterial 
growth experiments, each data point represents the mean value from at least thirty leaves. Similar results were obtained in two independent 
experiments. The bars represent standard deviation. (Student's t-test, * p<0,05). 
doi:1 0.1 371 /journal.pone.0099343.g002 



despite a significant basal activation in the absence of pathogen, 
ltf2 mutant plants were able to further activate SA-related defense 
genes in response to pathogen attack, and the response of the lif2 
mutant to P. syringae infection differed from the wild type only in 
the level of defense gene expression. Together, these results 
indicate that LIF2 is not necessary for the activation of SA- and 
JA-related defense genes. 

Response of the Iif2 sid2 double mutant to the biotrophic 
P. syringae pathogen 

Since neither the high level of PRl expression in the lif2 mutant 
nor the reduced susceptibility of lif2 to P. syringae infection were 
correlated with an increase in free SA content, we hypothesized 
that LIF2 may act downstream of SA production. To test this 
hypothesis, we crossed the lif2 mutant with the salicylic acid induction 
deficient2 (sid2/icsl) mutant. 

The rosettes of sid2 lif2-l double mutant plants were smaller 
than those of the two parental lines, sid2-2 and lif2-l, in both 
short-day (SD) and long-day (LD) conditions (Fig. 4A and S2). 
Interestingly, in SD conditions, similarly to lif2-l plants, the 
double mutant was early flowering (Fig. S2) and produced 
indeterminate ovaries (IDO) with an ectopic inflorescence, 
suggesting that the lif2-l mutation is epistatic to sid2-2. 

We then investigated the response of the lif2-l sid2-2 double 
mutant to pathogens. Bacterial growth of the DC3000 strain was 
reduced at 24 hpi (by about 10-fold) in the lif2-l sid2-2 rosette 
leaves compared with those of the wild type (Fig. 4B), similarly to 
lif2-l. These results suggest that the reduced susceptibility of lif2 is 
independent of SA biosynthesis. Interestingly, the level of PRl 



expression in lif2-l sid2-2 was similar to that of the wild type 
(Fig. 4C), suggesting that the decreased susceptibility of lif2 may 
involve an &4-independent defense pathway. 

Response of the Iif2 mutant to a necrotrophic fungal 
pathogen 

Since responses to pathogens are usually largely dependent on 
the pathogen lifestyle, we then assessed the defense response of lif2 
to the necrotrophic fungal pathogen S. sclerotiorum (strain S55). We 
used two A. thaliana accessions as controls, Rubezhnoe-1 (Rbz-1) 
and Shahdara (Sha), which were previously shown to be resistant 
and susceptible to S. sclerotiorum, respectively [29]. After inocula- 
tion, the symptoms were stronger in lif2 than in wild-type plants 
(Fig. 5A and B), revealing that lif2 is susceptible to S. sclerotiorum. 
This observation is consistent with the constitutive downregulation 
of JA-related genes in lif2 reported in this study, and the finding 
that the JA pathway is essential for resistance to this pathogen 
[29,30,31]. 

We then inoculated these plants with a second necrotrophic 
fungal pathogen, B. cinerea. Unexpectedly, we observed that ltf2 
was less susceptible to two virulent strains of B. cinerea, the wild- 
type virulent B05 1 0 and Bd90 strains, the latter of which was less 
virulent (Fig. 6A-C). These data suggest that the immune response 
of lif2 is independent of the pathogen's lifestyle. 

To decipher the molecular basis for the resistance of lif2 to B. 
cinerea, we assessed the expression of a number of regulators known 
to be involved in disease resistance to B. cinerea (Table 2). For 
instance, mutations in the AflTJ^ff transcription factor gene [32], 
ATGXRS13 [33], or the basic helix-loop-helix MYC2 transcription 
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Figure 3. Expression of genes involved in tKie JA- and SA-dependent signaling pathways in response to P. syringae inoculation. 

Seven-week-old rosette leaves of wild-type (WT) and lifl-l mutant plants were inoculated with buffer (mock) or with the bacterial pathogen DC3000. 
Leaves were collected at 24 hpi and gene expression was assessed by quantitative RT-PCR. The expression of genes involved in (A) JA- and (B) SA- 
dependent signaling pathways is shown. WT mock (white), WT DC3000 (light grey), Iif2-1 mock (grey), and Ilf2-1 DC3000 (black). 
doi:1 0.1 371/journal.pone.0099343.g003 



factor [34] led to increased resistance to B. cinerea. In our 
transcriptome data, ATGXRS13 was upregulated in lif2, while 
MTB46 was not affected, suggesting that these two genes are 
probably not involved. However, we found that MTC2 was 
downregulated in lif2 (Fig. 6D), which is in agreement with the 
decrea.sed susceptibility of lif2 to B. cinerea. We also evaluated the 
expression profiles of a few other candidate genes whose 
downregulation led to an increased susceptibility to B. cinerea 
(i.e., PAD2, CYP71B15/PAD5, BIKl, WRKTIS, WRKTSS, 
WRK40, and WRKYTO) [35,36,37,38,39,40]. Interestingly, only 
the expression of WRKT18 and WRKT33 was significantiy different 
in the two mutant alleles compared with wild-t>'pe or lif2-c plants. 
The loss-of-function of CATALASE 2 {CAT2) triggers pathogen 
defense responses and resistance [41]. Furthermore, the inhibition 
of NADPH oxidase and of other flavoprotein enzymes involved in 
oxidative stress limits B. cinerea infection [42] . We found that CA T2 
was downregulated in lif2, which is in agreement with the reduced 
susceptibility of this mutant to B. cinerea. Thus, CAT2 and the 
MTC2, WRKT18, and WRK33 transcription factors are interesting 
candidate proteins that participate in the lif2 response to B. cinerea. 

Glucosinolate content is altered in the Iif2 mutant 

Due to the observed enrichment in the GO term "glycosinolate 
metabolic process" (GO:0019757) in the lif2-l transcriptome 
(Table 1), and to the established role of glucosinolates (GLSs) and 
their breakdown products in plant defense [43] , we extracted and 
quantified twenty-one of these secondary metabolites (fi-thiogluco- 
side-N-sulfated oximes), belonging to three main GLS families 



(aliphatic, indolic, and benzoate GLSs) in mutant and wild-type 
seedlings (Table SI). The GLSs were analyzed by negative 
electrospray ionization liquid chromatography coupled with mass 
spectrometry (ESI-HPLC-MS). Globally, the GLS content was 
lower in the liJ2 mutant than in wild-type plantlets, with decreases 
observed in fourteen GLSs belonging to all three main GLS 
families (Fig. 7A). The levels of indol-3-ylmethyl glucosinolate 
(ISM), 4-methylsulfmylbutyl glucosinolate (4MSOB), 4- 
methylthiobutyl glucosinolate (4MTB), and 5-methylthiopentyl 
glucosinolate (.^MTP) were significantly decreased in the lif2 
mutant. Interestingly, the indolic GLS family was less affected by 
the lif2 mutation, with only a decrease in indol-3-ylmethyl 
glucosinolate (glucobrassicin, ISM) being observed. As a control, 
we quantified the GLS content in the Ihpl mutant (Fig. 7B), which 
is known to have an altered glucosinolate profile [15]. In Ihpl, 
significantly increased levels were obscr\-cd for six GLSs (i.e., 
8MSOO, SMTP, 4MSOB, 5MSOP, 6MSOH, and 7MSOH), 
whereas only one GLS, the indolic 4MOI3M GLS, had decreased 
levels relative to the wild type. Similar results were obtained using 
Ihpl seeds [15]. Next, we quantified the levels of raphanusamic 
acid (RA), an important breakdown product of GLS that forms 
during the biotic stress response [44]. PEN2 myrosinase catalyzes 
the formation of RA from the ISM or 4MOI3M substrates. We 
observed a significant reduction in RA in both lif2 (0.91 ng/mg, 
n = 7-8, p<0.01) and Ihpl (0.89 ng/mg, n = 7-8, p<0.001) 
compared with wild-type plants (1.39 ng/mg). 

Therefore, lif2 has a GLS profile that is distinct from that of 
wild-type and Ihpl seedlings. Thus, the control of GLS metabolism 
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Figure 4. The response of the Iif2-1 sid2-2 double mutant to P. syringae inoculation. (A) Seven-week-old plants grown In short-day 
conditions. Scale bar, 1 cm. (B) DC3000 bacterial growth at 24 hpi. Each inoculated leaf was ground in IVIgCl2 and the bacterial suspension was then 
diluted and plated on solid medium. Twenty-four leaves were analysed per genotype. Stars indicate a significant difference from wild-type (WT) 
plants (Mann and Whitney, * p<0.05). Experiments were repeated twice and gave similar results. (C) The relative expression of PR1 in untreated 
rosette leaves of seven-week-old plants. 
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seems to be dependent on LIF2 and partially independent of the 
LHP1/LIF2 interaction [12]. Alterations to the GLS pathways 
may also contribute to the responses of the lif2 mutant to 
pathogens. 

Molecular signature of the lifl mutant 

In an attempt to identify the key genes involved in the lif2 
response to biotic stress, we further analyzed the expression of 
transcription factors (TFs), due to their prominent role in 
regulating gene expression. Among the 66 TFs deregulated in 
the lif2-l transcriptome [12], we noticed that the APETALA2/ 
ethylene-responsive element-binding protein (AP2/EREBP), 
WRKY, MYB, and basic helix-loop-helLx (bHLH) transcription 
factor families were overrepresented (Table S2) (Fig. 8A). We also 
noticed a strong bias towards downregulation of genes of these 
families (Fig. IC). For example, 11 of the 13 deregulated genes of 
the AP2-EREBP family were downregulated in lif2 (Fig. 8). Out of 
the four subfamilies of AP2/EREBP (i.e., AP2, RAV (related to 
ABI3/VP1), dehydration-responsive element-binding protein 
(DREB), and ERF) [45,46], only the two largest ones, the DREB 
and ERF subfamilies, were deregulated in the lif2-l mutant. The 
plant-specific WRKY transcription factors are key regxilators of 
stress and plant immune responses [47], whereas the NAG TFs are 
involved in both development and the abiotic and/ or biotic stress 
responses. Two stress-responsive NACs (SNACs) that were 



recently described [48,49] were deregulated in the lif2-l 
transcriptome. In the MYB and bHLH families, MYC2, which 
encodes a key defense transcription factor involved in JA 
responses, was deregulated, (Fig. 6D) as well as the JA- 
ASSOCIATED MrC2 UKEl gene (JAMl). 

Furthermore, we noticed that 6 genes associated with the GO 
term "negative regulators of defense response" (GO:0031348) 
were downregulated in the lif2-l transcriptome. These genes 
might participate in the liJ2 stress response (Fig. 8B). Among these 
genes, JARl represses the SA-dependent signaling pathway; 
PEN3, which is an ABG transporter, restricts pathogen prolifer- 
ation in the hosts; and BAPl is a negative regulator of 
programmed cell death and is involved in membrane trafficking 
in response to external conditions [50]. These data suggest that 
LIF2 affects a variety of defense-related pathways. 

The 38 TFs that were both deregulated in the lif2-l rosette 
transcriptome (Table S2) and belong to the AP2/EREBP, WRKY, 
NAG, MYB, or bHLH family, as well as the 6 negative regulators 
(Fig. 8B) represent a sort of molecular signature of the lif2 mutant 
phenotype. We thus wondered how these 44 genes were 
deregulated in wild-type genetic backgrounds in response to biotic 
stress. To address this question, we selected 62 transcriptome 
profiles from plants treated with pathogens compared with those 
that were not, in the Genevestigator database [51] and extracted 
the expression of these 44 genes. We then performed a hierarchical 
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Figure 5. LIF2 is involved in the plant defense response to 5. sclerotiorum. Leaves of four-week-old plants were Inoculated with S. 
sclerotiorum strain S55. The Iif2 alleles and the complemented Iif2-c line in the Col-0 background were analysed. A. thaliana Col-0, Rubezhnoe-1 (Rbz- 
1) (more resistant than Col-0), and Shahdara (Sha) (more susceptible than Col-0) accessions were used as controls. (A) Symptoms at 7 dpi. (B) The 
disease score was evaluated for each line at 7 dpi. Means and standard deviations were based on at least twenty plants per line. A significant 
difference in susceptibility relative to that of the Col-0 accession is indicated with an asterisk (Kruskal and Wallis's test, * p<0.05). 
doi:1 0.1 371 /journal.pone.0099343.g005 



clustering analysis using of the MultiExperiment Viewer applica- 
tion [52], and identified three main gene clusters, which are 
characterized by specific co-expression profiles of MTC and MTB 
(Tig. 8C). The small Cluster I (containing MrC4/MrB28/MrB48) 
contained genes that are upregulated in lif2-l and downregulated 
in response to most pathogen treatments. The large Cluster II 
(containing MYC2/MYB15) contained genes downregulated in 
lifl-l and up- or downregulated in response to pathogen 
treatments, with no obvious link to the biology of the pathogens. 
Finally, Cluster III (containing MTC7E/MTB34) mainly contained 
genes repressed in lif2-l and in response to most biotic stresses, 
except in response to P. syringae, and thus, which are probably 
associated with the JA pathway due to the presence of JAR 1. In 
conclusion, our clustering analysis revealed that (i) the expression 
profiles of the 44 selected genes involved in the response to 
pathogens were diverse, (ii) the MYC and MYB genes were co- 
expressed in different combinations, and (iii) that lif2 exhibited a 
unique expression pattern, with opposite profiles to most of the 62 
analyzed transcriptome profiles, which might be in line with its 
specific response to pathogens. 

Discussion 

Plants have evolved both constitutive and induced defense 
mechanisms to counteract pathogen attacks. Most of the induced 
mechanisms are based on the production of a complex repertoire 
of plant metabolic compounds and signaling hormones, which 
activate appropriate defense pathways, based on the nature of the 
pathogens [53,54]. Salicylic acid, a phenolic phytohormone, and 
the jasmonic acid phytohormone, a polyunsaturated fatty acid 
derived from ot-linolenic acid, are two key molecules in plant 



defense mechanisms [55]. SA, which also possesses medicinal 
properties [24] and plays a central role in animal immunity [56], 
and its derivatives, have crucial roles in defense against biotrophic 
pathogens [57]. They also regulate various aspects of abiotic stress 
responses, plant growth, and development. The JA and SA defense 
signaling pathways have been shown to cross-communicate and to 
be mostly antagonistic, providing plants with a regulatory potential 
to fine-tune their defense reaction depending on the type of 
pathogen encountered [22,23]. SA and JA are mainly involved in 
the response to biotrophic and necrotrophic pathogens, respec- 
tively. However, complex crosstalk exists between the hormone 
signaling pathways involved in plant responses to pathogens [26]. 
Secondary metabolites, such as glucosinolates, are also involved in 
the biotic stress response in plants [43]. 

Only a few RBPs have emerged as regulators of the biotic stress 
response in plants [8,58,59,60]. For instance, AtRBP-DRl, a 
putative component of the resistance protein RPS2 complex, 
positively regulates defense responses mediated by SA in A. thaliana 
[61]. The conserved MOS2 protein, a putative RBP with both G- 
patch and KWO domains, is required for resistance against P. 
syringae pv. maculicola ES4326 and Peronospom parasitica Emoy2 [62], 
whereas overexpression of the Gossypium hirsutum GhZFPl protein, 
which contains a CCCH-type zinc finger present in some RBPs, 
has been shown to enhance both disease resistance to the 
Rhizoctonia solani fungus and salt stress tolerance [63]. Over- 
expression of the Capsicum annuum RNA-binding protein 1, 
CaRBPl, in^. thaliana conferred reduced susceptibility to infection 
by the biotrophic oomycete Hyaloperonospora arabidopsidis [9]. 
Except for the hnRNP AtGRP7 [9], the RBP-mediated defense 
mechanisms remain poorly understood. 
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Figure 6. LIF2 is involved in the plant's susceptibility to B. cinerea. (A-C) Six-weel<-old plants were infected with a mycelium plug of the 
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doi:1 0.1 371 /journal.pone.0099343.g006 



Here, we studied the role of ^. thaliana RBP LIF2 in the plant's 
response to biotic stress. We demonstrated that LIF2 loss-of- 
function results in pathogen resistance, in a manner that is 
independent of the pathogen's lifestyle or infection strategy. 
Indeed, the lif2 mutant was less susceptible to both the hemi- 
biotrophic bacterial P. syringae pv. tomato pathogen and the 
necrotrophic ascomycete B. cinerea, two pathogens with substantial 
impacts on agronomical production due to their wide range of host 
specificities [17,18]. However, lif2 was more susceptible to the 
necrotrophic ascomycete S. sclerotiorum, which attacks more than 
400 plant species around the world. 

The three pathogens used in this study induce different plant 
defense pathways and their associated defense molecules: the JA 
and ethylene pathways, and possibly the SA pathway, are 
activated in response to S. sclerotiorum infection [29,30]; the SA 
pathway is key in the plant's response to P. syringae infection; and 
JA and camalexin are involved in the response to B. cinerea 



infection [64] . Our analysis of the lif2 mutant revealed that SA- 
related genes [ICSl / SID2, PRl) were constitutively activated. The 
upregulation of PRl was in agreement with the decreased 
susceptibility of lif2 to P. syringae, and this response was also 
observed in other mutants with reduced susceptibility to P. syringae, 
such as cprl and cpr5 [65]. However, the constitutive activation of 
the SA biosynthesis pathway did not lead to an overaccumulation 
of free SA in the lif2 mutant. SA undergoes various chemical 
modifications, which alfect its activity, catabolism, transport, and 
storage [66] . The accumulation of free SA is fine-tuned by the 
transcriptional regulation of genes involved in SA biosynthesis, but 
also by the modulation of enzymes that modify SA [24] . Our data 
suggest that genes of the SA pathway may undergo post- 
transcriptional regulation or that a homeostasis mechanism may 
control free SA accumulation and maintain it at wild-type levels in 
lif2. Furthermore, the SA production pathway is not fuUy 
elucidated in A. thaliana, since no isochorismate pyruvate lyase, 
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Table 2. Genes involved in susceptibility to 6. dnerea. 



Name 


AGI 


Mutant phenotype 


Reference 


GRXS13 


AT1G03850 


Decreased susceptibility 


La Camera et al. 201 1 


MYC2 


AT1G32640 


Decreased susceptibility 


Lorenzo et al. 2004 


MY846 


AT5G 12870 


Decreased susceptibility 


Ramirez et al. 201 1 


BIKl 


AT2G39660 


Increased susceptibility 


Veronese et al. 2006 


con 


AT2G39940 


Increased susceptibility 


Thomma et al. 1999 


JARl 


AT2G46370 


Increased susceptibility 


Thomma et al. 1999 


PAD2 


AT4G23100 


Increased susceptibility 


Ferrari et al. 2003 


PADS 


AT3G26830 


Increased susceptibility 


Ferrari et al. 2003 


WRKY18 


AT4G31800 


Increased susceptibility* 


Xu et al. 2006 


WRKY40 


ATI G80840 


Increased susceptibility* 


Xu et al. 2006 


WRKY33 


AT2G38470 


Increased susceptibility 


Zheng et al. 2006 


WRKY70 


AT3G56400 


Increased susceptibility 


AbuQamar et al. 2006 



doi:1 0.1 371 /journal.pone.0099343.t002 



which is required to convert isochorismate into SA in bacteria, has 
been identified to date. A second SA biosynthesis pathway, which 
is regulated by phenylalanine ammonia lyase (PAL), has been 
proposed to exist in plants [57,67] . In lif2, the discrepancy between 
the constitutive upregulation of the SA/ICSl -dependent pathway 
and the lack of SA accumulation may also suggest some crosstalk 
between the ICSl and PAL pathways to fme-tune SA production. 
Finally, the ltf2 sid2 double mutant was less susceptible than wild- 
type plants to P. syringae, despite exhibiting a PRl expression level 
similar to that of wild-type plants. Therefore, the reduced 
susceptibility of the lif2 mutant to P. syringae was independent of 
SA in non-challenged plants, suggesting that LIF2 participates in 
plant defense via a novel defense pathway that is independent of 
the SA-signaling defense pathway. 

Furthermore, we showed that the GLS profile of the lif2 mutant 
was altered. These changes might indirectly participate in the 
increased resistance phenotype of lif2. Indeed, GLSs function in 
the plant's defense response to herbivores and fungal pathogens 
[44,68,69,70] . Ward et al. (20 10) suggested that GLSs do not play a 
direct role in the response to bacteria, but possibly have an 
intermediary role by influencing some defense pathways. Indeed, 
after infection of ^. thaliana with P. syringae DC3000, a significant 
reduction in ISM and 4MOI3M was observed [71]. In the lif2 
mutant, a low level of ISM was observed, but no change was 
observed in 4MOI3M in the absence of bacterial challenge. The 
global decrease in GLS and RA contents observed in lif2 might 
thus contribute to the control of the basal defense pathways. 
Interestingly, JA is a regulator of GLS gene expression and of GLS 
accumulation via multiple pathways [72,73]. 

Several TFs have recendy been shown to be involved in GLS 
gene regulation, including MrB28, MrB29, and MrB76, which 
are key regulators of the aliphatic-GLSs, and MYB34 and A4TB51, 
which are regulators of the indole-GLSs. Furthermore, MYC2, 
MYC3, and MYC4 can form protein complexes with all known 
GS-related MYBs to regulate GLS biosynthesis [74]. In the lif2 
mutant, MYB28 is downregulated, whereas MYB34 is upregulated, 
and the levels of members of the two GLS families were decreased. 
However, MYC2 is downregulated in the lif2 mutant. Further- 
more, several TFs related to defense responses are deregulated in 
lif2. WRKY33 is an essential transcription factor in the defense 
against B. cinerea that acts by controlling the expression of genes 
involved in redox homeostasis, SA signaling, and camalexin 



biosynthesis and thus affecting the SA-JA balance [40]. The 
activation of WKRY33 observed in lif2 may have downstream 
effects, which may play a role in the resistance of lif2 to B. cinerea. 
Therefore, it is likely that multiple components of different defense 
pathways contribute to the primed state of lif2 in the absence of 
pathogen, and to its reduced pathogen susceptibility. Interestingly, 
LIF2 interacts with the chromatin-associated protein LHPl, a 
subunit of the Polycomb Repressive Complex, which interacts 
with numerous genomic sites and regulates their expression 
[12,14]. Some defense-related genes are present among the 
LHPl targets. Whether LIF2 acts coordinately with LHPl at 
these defense-related loci to control them constitutes an interesting 
research question, as the role of chromatin proteins in plant 
immunity is poorly documented [75,76,77,78,79]. For instance, 
PIE 1 , a member of the SWR 1 subfamily, was shown to negatively 
regulate plant defense [80], whereas loss-of function of the 
SIRTUIN2 histone deacetylase, a homolog of the yeast Silent 
information repressor 2 (Sir2) protein [81] and of SDG8 histone 
methyltransferase [82] alter plant-pathogen responses. It was 
proposed that chromatin modifications may participate in defense 
priming in plants [83]. Thus, our study highlights an emerging role 
for the LIF2 chromatin-associated protein in biotic stress 
responses, with a putative suppressor function in plant immunity 
(Fig. 9). 

Finally, LIF2 is also involved in the maintenance of growth and 
cell determination during floral development [12]. Therefore, our 
study illustrates the trade-off existing between plant development 
and plant defense against various enemies in a changing 
environment [84,85] . LIF2 may thus regulate the balance between 
development and plant immunity by limiting the cost associated 
with the plant defense response (Fig. 9). It was suggested that 
priming of SA-related defense responses greatiy enhances disease 
resistance and plant fitness, but diminishes fitness in the absence of 
pathogens [85] . Furthermore, since both developmental and plant 
defense processes require adaptation to environmental conditions, 
sharing common elements with dual functions may allow for better 
fine-tuning [86]. Since the human linRNP-Q^ homolog of LIF2, 
SNAPl /SYNCRIP, also acts as a cellular RBP and as a suppressor 
of human immunity against virus infection, our data highlight a 
conserved role for hnRNP-Q_ in eukaryote immunity. The dual 
function of hnRNP-Q, proteins in development and defense is 
likely to involve key conserved molecular events in eukaryotic cells. 
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Figure 7. Glucosinolate (GLS) contents in li'f2-7 and Ihp1-1 young seedlings. The GLSs were quantified In tifl-t (A) and /hpI-7 (B) using ESI- 
HPLC-MS. The full names of the GLSs are listed in Table SI. 
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Figure 8. Stress-response genes deregulated in the //^transcriptome. (A) The distribution of dereguldte(d TPs in Iif2. (B) Deregulated genes 
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data. (C) Hierarchical clustering analysis performed using the MultiExperiment Viewer application, the Pearson correlation as current metric, and 
complete linkage clustering as the linkage method. The expression profiles of 38 deregulated TPs and 6 negative regulators in different biotic 
conditions were used in this analysis. Three gene clusters (I to III) were identified. 
doi:1 0.1 371/journal.pone.0099343.g008 
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Unknovm mechanisms 



Figure 9. IVIodei for tKie functions of the hnRNP-Q LiF2 protein. 

LIF2 may regulate the balance between development and plant 
Immunity by minimizing the energy cost of plant defense. 
doi:l 0.1 371/journal.pone.0099343.g009 

Identifying the. associated RNA partners of LIF2 RBP may 
elucidate the underlying mechanism. 

Materials and Methods 

Plant materials 

Amhidopsis thaliana lines used in this study were in the Columbia 
(Col-O) accession. The homozygous lif2-\ and lif2-3 mutants were 
described in [12]. The lif2-c line is a lif2-l mutant complemented 
with the LIF2 genomic region under the control of a 3-kb 
promoter region. The sid2-2 [87] mutant was crossed with the lif2- 
1 allele. Double mutants were selected by PGR using specific 
primers (Table S3, [12]). 

For P. syringae, B. cinerea, and S. sclerotiorum pathogen response 
analyses, plants were grown on soil in growth chambers, under 
controlled short-day conditions (8 hours light/ 16 hours dark), at 
20°C and 70% hygrometry. 

Primers 

All primers are described in Table S3. 
Bacterial inoculation 

The virulent strain P. syringae pv. tomato DC3000 [Pst DC30()0) 
[88] and the avirulent strain P. syringae pv. tomato DC3000, 
carrying the avirulent gene avrRpml {Pst DC 3000 avrRpml) [28], 
were grown at 28°C on King's B solid medium plates containing 
the appropriate antibiotics. Briefly, a fresh bacterial suspension 
was scraped off the growth plate, and resuspended in 10 mM 
MgClj at an OD600 of 0.01 (at a concentration of 10^ cfu ml"'). 
The bacterial suspension was infiltrated into the abaxial side of 24 
to 30 rosette leaves of seven-week-old plants (6 to 10 plants per 
genotype) using a 1-ml syringe without a needle. Control 
inoculations (Mock) were performed using 10 mM MgCl2. 
Bacterial growth was assessed 24 hours after infiltration. Whole 
leaves were homogenized in 500 \A of 10 mM MgCl2. After serial 
dilutions, homogenates were plated on King's B medium and 
incubated at 28°C for 2 days before colonies were counted. The 
surface area of each leaf was measured using ImageJ software in 
order to calculate cfu/cm^. Day 0 titers ranged from logio - 3.6- 
3.9. Leaves were photographed 5 days post-infiltration (dpi). All 
experiments were repeated at least twice with similar results. 

Fungal pathogen inoculation procedure 

Four-week-old plants were inoculated with discs carrying S. 

sclerotiorum (S55 strain) mycelia, as previously described [29]. Plant 
inoculations were performed in a Helios 1200 Phytotron under a 
light period of 9 h at 22°C and 80% relative humidity during the 



first 3 days and 66% thereafter. At least 20 plants of each genotype 
were inoculated. The disease score of each line was evaluated 
seven days post-inoculation. Two independent experiments were 
performed with similar results. 

Leaves of seven-week-old plants were inoculated with mycelial 
plugs of the wild-type strain B. cinerea B0.5 1 0 and Bd90 (diameter of 
3 mm), as previously described [89] . Infected leaves were placed in 
a Petri dish under high humidity and incubated at 21°C. Lesion 
surfaces were measured daily for three days using ImageJ software. 

Quantification of hormonal contents 

Rosettes of seven-week-old plants were frozen immediately after 
harvest and ground in liquid nitrogen. Four pools of 6 rosettes 
were collected. Frozen material (100 mg) was extracted with 3 ml 
of acetone/water/acetic acid (80/19/1, v:v:v) containing the 
following stable isotope labelled internal standards: 10 ng [4-2H] 
ABA (NRC-CNRC Plant Biotechnology Institute, Saskatoon, 
Canada), 50 ng [4-2H] sahcylic acid (Olchemlm, Olomouc, Czech 
Republic), 1 ng [5-2H] jasmonic acid (CDN Isotopes GIL 
Cluzeau, Sainte Foy la Grande, France), and 10 ng [6-13G] 
indole-3-acetic acid (Cambridge Isotope Laboratory, Andover, 
MA). The extract was vigorously shaken for 30 s, sonicated for 
1 min at 25 Hz, shaken for 10 min at room temperature, and then 
centrifuged (8230 g, 4°C, 15 min). The supematants were 
collected and the pellets were extracted again with 1 ml of the 
same extraction solution, and then vigorously shaken (1 min) and 
sonicated (1 min, 25 Hz). Following centrifugation, the two 
supernatants were pooled and dried. The dry extract was dissolved 
in 140 |i,l acetonitrUe/water (50/50, v/v), filtered, and submitted 
to analysis by HPLC-electrospray ionisation-MS/MS (HPLG-ESI- 
MS/MS). The compounds were introduced into the ESI source 
using a Waters 2695 separation module (Alliance; Waters, Milford, 
MA, USA) equipped with a Waters 2487 dual UV detector. 
Separation was achieved on a reverse-phase column (Uptisphere 
C18 ODB, 150*2.1 mm, Interi him, Montlugon, France), using a 
flow rate of 0.15 ml min ' and a binary gradient as follows: (A) 
acetic acid 0.1% (v/v) and (B) acetonitrile. The solvent gradient 
was programmed as follows: 0-5 min, 20% A; 5-15 min, 50% A; 
15-30 min, 100% B; and 30-42 min, 20%. The analyses were 
performed on a ^\'aters Quattro LC Triple Qiiadripole Mass 
Spectrometer (Waters) operating in a Multiple Reaction Moni- 
toring (MRM) scanning mode. The instrumental parameters were 
set as follows: capillary, 2.70 kV (negative mode); extractor, 3 V; 
and source block and desolvation gas temperatures, 120°C and 
350°G, respectively. Nitrogen was used for the nebulization and 
desolvation (77 L h ' and 365 L h respectively), and argon was 
used as the collision gas at 2.83 10 ^ mbar. For a 5-nL injection 
volume of sample prepared and reconstituted in 140 nl of 50/50 
acetonitrile/H20 (v/v), the limit of detection (LOD) and limit of 
quantification (LOQ) were extrapolated for each hormone from a 
calibration curve and sample using the Quantify module of 
MassLynx (version 4.1 software). The parameters used for MRM 
quantification and the LOD and LOQ are listed in Tables S4 and 
S5, respectively. The amount of JA was expressed as a ratio of 
peak areas (209>62/214>62) per fresh weight, due to impurities 
contained in the D5-JA standard. 

Glucosinolate extraction and quantification 

Glucosinolate analyses were performed on 15-day-old in vitro 
seedlings grown under long-day conditions. The plant material 
was rapidly collected, frozen, and lyophUized. LyophUized 
material (30-50 mg) was ground in liquid nitrogen. Then, 
0.2 mg of sinigrin hydrate standard (Sigma Aldrich Ref 85440) 
was added as an internal tracer for recovery and analytical 
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purposes. The samples were extracted using 2.5 ml of extraction 
solution (75% acetonitrile/25% water solution), vigorously ho- 
mogenized for 5 min in a Polytron homogenizer (Fischer 
Scientific), and centrifuged (8000 g, 20°C, 15 min). The superna- 
tants were collected and the pellets were extracted with 1 ml of the 
extraction solution and sonicated (15 min, 25 Hz). Following 
centrifugation, the supernatant was pooled with the first superna- 
tant. The extracts were then evaporated to dryness using a 
Thermo Savant Speedvac overnight, at ambient temperature. A 
water/acetonitrile solution (95/5) was added to the dried extracts, 
whi[:h were filtered and diluted to the third with water/ acetonitrile 
(95/5), prior to negative electrospray LC-MS analysis on an 
Alliance 2695 system coupled to Quattro LC (Waters). Glucosi- 
nolates were identified by retention time, mass, isotopic pattern, 
and fragment ions [90,91]. The concentrations of the metabolites 
of interest were quantified using the sinigrin response, as 
previously described [92]. Three independent extraction analyses 
were carried out per biological experiment and two to three 
biological replicates were performed. Biochemical information 
about the different classes of GLSs and their biosynthesis pathways 
can be found in various databases (e.g., http:/ /www.genome.jp/ 
kegg-bin/show_pathway?ko00966-l-C08417 or the AraCyc web- 
site). 

Gene expression analysis 

Total RNA was isolated from various tissues, using the RNeasy 
Plant Mini Kit (QIAGEN), and treated with RNase-free DNasel 

(Invitrogen). Reverse transcription (RT) reactions were performed 
with Superscript II reverse transcriptase (Invitrogen), according to 
the manufacturer's instructions. Quantitative real-time PGR 
(qPCR) was performed on Eppendorf Mastercycler® ep realplex 
(Eppendori) using MESA FAST qPCR MasterMix Plus for 
SYBR® Assay (Eurogentec), as per the manufacturer's instructions. 

Supporting Information 

Figure SI LIF2 expression in response to P. syringae 
DG3000 and DC3000 avrRptnl inoculations, hpi, hour 
post-infection. 
(PPTX) 

Figure S2 Phenotypes of the liJ2-l sid2-2 double mutant. 

(A) Rosettes of 46-day-old plants grown in long-day (LD) 
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